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Abstract. The fluctuations of spacetime geometries at finite temperature are 
evaluated within the linearized theory of gravity. These fluctuations are described 
by the probability distribution of various configurations of the gravitational field. 
The field configurations are described by the linearized Riemann-Weyl tensor 
without any reference to the metric. The probability distribution of various con¬ 
figurations is described by the Wigner functional of the gravitational field. It 
has a foam-like structure; dominant configurations are those with large changes 
of geometry at nearby points. In the high-temperature limit one obtains the 
Bolzmann distribution that enables one to identify the expression for the total 
energy of the gravitational field. The appearance of the same expression for the 
total energy when the gravitational field is treated as a collection of gravitons and 
as the high-temperature limit of the Wigner functional proves the consistency of 
the whole procedure. Striking differences are found between the fluctuations of 
the electromagnetic field and the gravitational field; among them is the divergence 
in the gravitational case of the probability distribution at zero temperature. This 
divergence is of the “infrared type” because it occurs in integrals over the wave 
vector at small k. 
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The main goal of this paper is to describe the quantum fluctuations of the gravitational 
field at finite temperature in the linearized gravity. This is done with the use of the 
Wigner function; a tool that has never been used in this context before. In order to 
reach this goal I shall introduce the quantization procedure that has the advantage 
of involving only the geometrical, gauge invariant elements. This is accomplished in 
the framework of the spinorial formulation providing another evidence for the special 
role played by spinors and justifying once again the praise of spinors expressed by 
Andrzej Trautman [T]: “One is tempted to say that the world around us, and life in 
particular, are so rich because Nature found it convenient to use, among its building 
blocks, entities requiring spinors in their description.” 

The original Wigner function is a function of positions and momenta of quantum 
particles. However, there is natural generalization of this concept. Namely, we 
may replace the canonically conjugate particle variables by their field-theoretic 
counterparts. This was done for the scalar field in [5] and for the electromagnetic field 
in [3] . Upon this generalization, the Wigner function becomes the Wigner functional 
whose arguments are the field configurations. 

Quantum properties of fields manifest themselves, in particular, in the field 
fluctuations present even in the vacuum state. In quantum electrodynamics these 
fluctuations lead to observable effects (Welton explanation of the Lamb shift, Casimir 
effect, photon shot noise). In classical physics statistical field fluctuations are described 
by the probability distribution of various field configurations. In quantum physics this 
simple description fails, due to the uncertainty relations. However, when the Wigner 
function is positive it may serve as a very good substitute for the standard classical 
distribution function. 

Statistical properties of the vacuum fluctuations of the gravitational field are 
embodied in the probability distribution that assigns relative weights to different 
geometries. Exact solution of this problem is a hopeless task since it would require 
full-fledged quantum theory of gravity. However, an approximate solution can be 
obtained in linearized gravity—after all, our Universe is mostly flat. In this approach 
I assume that the linearized gravitational field can be quantized just like any other free 
field. Every free field can be viewed as a collection of uncoupled harmonic oscillators. 
Since the Wigner function for the thermal state of harmonic oscillators is Gaussian, 
it can serve as a bona fide probability distribution. 

I do realize that linearization is a crude and questionable procedure characterized 
long time ago by Roger Penrose in highly critical terms [1] “... if we remove the life 
from Einstein’s beautiful theory by steam-rollering it first to flatness and linearity, then 
we shall learn nothing from attempting to wave the magic wand of quantum theory over 
the resulting corpse”. Nevertheless, I did make full use of an important feature of the 
Penrose approach—the correspondence (fully underscored in the spinorial formulation) 
between the electromagnetic field and the Riemann tensor. 

In a recent paper Freeman Dyson questioned altogether the widely accepted 
view that the gravitational field is just another field waiting to be quantized. He 
raised the possibility that “the gravitational field is a statistical concept like entropy 
or temperature, only defined for gravitational effects of matter in bulk and not for 
effects of individual elementary particles”. If this indeed would be so and the gravitons 
would not exist, my analysis will loose its foundation. However, if you cannot do what 
you would like, you should like what you can do. 
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The analysis of the gravitational Wigner functional is greatly simplified if we 
take full advantage of the analogy between the electromagnetic field tensor and 
the Riemann tensor R^yXp in linearized gravity. Despite this close analogy, however, 
the differences between the electromagnetic and gravitational cases are substantial. 
The most striking difference is the divergence of the gravitational Wigner functional 
in the limit of zero temperature which might signify some problems with the ground 
state. I will show that the source of these problems is found in the singular behavior 
of Riemann correlators. 

In order to make the analogy between the electromagnetic field and the 
gravitational field unambiguous we need an anchor point that will enable us to convert 
the loose analogy into precise mathematical formulas. I choose as this anchor point 
the universal formula for the energy valid for any free field expressed in terms of 
annihilation and creation operators. For all massless particles this formula can be 
written in the form: 



a\{k)ax{k), 


( 1 ) 


where k = |fe| and A = ± is the sign of helicity. I have chosen the relativistic 
normalization of the operators so that their commutation relations have the form: 

ax{k),a\,{k') 

This normalization leads to the relativistically invariant volume element in ([1]) with k 
in the denominator. 

Quantum statistical properties (Planck spectrum, Einstein A and B coefficients, 
etc.) are, of course, the same for the gas of photons and (linear) gravitons. Both 
particles are massless bosons with two helicity states. Large differences appear only 
in the description in terms of fields localized in space. 

The analogy between electromagnetism and gravity will be build on the universal 
formula (ED but the details will be most clearly seen in the spinorial formalism of 
relativity theory [B]. In this formalism only the true degrees of freedom will appear 
and there is no need for any gauge-fixing conditions. 

In the next section I briefly summarize the properties of this formalism needed 
for my analysis. In Sections [3] and S] I describe the quantized electromagnetic and 
the linearized gravitational field within the spinorial formalism. In Sections [5] and [B] 
I evaluate the Wigner functionals for the thermal states of the electromagnetic and 
gravitational fields and I discuss the physical properties of these states. In the last 
Section I discuss the problems with the ground state of the quantized gravitational 
field. In [Appendix A| I prove an important equation that relates the number of 
particles operators to the spinorial field operators. Finally, in [Appendix B| I use the 
Bel-Robinson tensor to fix the numerical coefficient in the normalization factor for the 
gravitational field. 


= 6xx'kS^^\k-k'). 


( 2 ) 


2. Spinorial formalism 

The most convenient representation of the fields describing massless particles is 
in terms of symmetric spinors. The connection between the spinors and the 
electromagnetic field tensor and the Riemann tensor will be given in Sec. [3| and 01 
respectively. Spinor indices will be denoted by capital letters and those for conjugate 
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spinors by dotted letters. We shall need the following 2x2 matrices: 


\9‘ 


mb'! 


AB 


AB 




AB 


= CAB, 


(3 a) 


- iT ^AB _ j 0 1 

J " - V -1 0 

5 ^ (3&) 

where tr’s are the Pauli matrices. My spinor conventions are those of Ref. m- They 
differ slightly from the conventions of Ref. [ 6 ] . The spinor indices take on the values 
( 0 , 1 ) and they are raised and lowered as follows: 

(t>^ = (j)B = (t>^CAB- (4) 

The spinorial wave function 4>ab...l{x) obeys the wave equation of the same form for 
all massless particles [5], 

g^^^d^(t)AB...L{x) = 0. (5) 

The number of indices is equal to 2H ; twice the absolute value of the helicity. 

The Fourier representation of the general solution of this equation, 

cUk 


/ d’^k 

(2^)3/2^ ^a«b ■■■bl [f+{k)e-^^'^ + f*_{k)e^^'^] 


( 6 ) 


expresses the decomposition of the field into harmonic oscillators. The presence of 
the volume element d^k/k underscores the relativistic content of this formula. The 
spinors ka are related to the integration variables k through the formulae: 

= k^, = \g>^^^k^. (7) 

The equations 0 do not determine the overall phase of ka- However, this phase is 
not significant because it can be absorbed by a change of the phases of the amplitudes 
f±{k) in classical theory or the phases of the annihilation and creation operators in 
quantum theory. A convenient choice of the spinor ka is: 

{ka} = —^ } {k,^-iky,k-k^}. ( 8 ) 

yj2{k - kz) 

The wave equation JS]) is satisfied due to the relations: 


yKA = g^'^^K^g^^ kbka = 2e^^ k 


B 


c^^ kakb = 0 . 


(9) 

The field equations retain the same form ([5]) also for the field operators 4’ab...l{x) 
in quantum field theory. In this case the amplitudes f±{k) in the expansion (jS]) into 
plane waves must be replaced by the annihilation and creation operators of particles 
with positive and negative helicity, 

d^k 


4>AB...Lix) = 7 


-kakb ■ ■ - Kl 


a+(fc)e-*'=-" + al(fc)e 


ik-: 


( 10 ) 


(27r)^/2fc 

The prefactor 7 is essential because the field operator ^ab...l{x) carries the 
dimensionality of the corresponding physical field while the dimensionality of the 
integral on the right-hand side is fixed by the canonical commutation relations ([ 2 ]). 
It follows from these relations that the annihilation and creation operators have the 
dimension of length. Therefore, the dimension of the integral is 1/length^^^. 

In what follows I shall need the following equality proved in [Appendix A[ 


^h{k) 


a{{k)ax{k) = 




OAAgOBB 


g (t>AB...L{r' 


( 11 ) 
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where the symbol f stands for T,\ f dfk/k, the function h{k) must be invariant under 
reflections, h{k) = h{—k), and 


h{r) 



ik-r Kk) 

ry2k‘iH-l 


( 12 ) 


For odd functions, h{k) = —h{—k), the contribution from negative helicity on the left 
hand side enters with the negative sign. In the formula (HU and also in (HU), (HZD, 
m, and (IA.3I) normal ordering of annihilation and creation operators is implied. 

The general formulation will now be applied to the electromagnetic field and 
then to the gravitational field. The corresponding annihilation and creation operators 
will be denoted by (c, c^) and {g,g^), respectively, while (a, a^) will be used in the 
general case. The well established electromagnetic case will serve as a guide in the 
construction in the gravitational case. I shall often use interchangeably the classical 
fields and their quantum counterparts but it should be clear from the context what is 
meant. 


3. Quantized Maxwell theory 


Classical electromagnetic field may be described by the second-rank symmetric spinor 
4>ab{x). The corresponding field operator (pAsix) is connected with the electro¬ 
magnetic field operator through the formulas (dotted indices for operators imply 
Hermitian conjugation): 

(13) 

(14) 


^ab{x) = 

U^ix) = ^ 


AS, 


4>ab{x) -\- <^ab 


The numerical coefficients in these formulas follows from the convention adopted in 
[5]. The value of the electromagnetic prefactor 'je appearing in the formula: 

dfk 


(pAB^x) = 'Ye 


-kabb 


c+(fc). 


-ik-: 


z'_{k)t 


(15) 


(27r)3/^fc 

can be found by comparing the expression (HD for the energy operator of the 
electromagnetic field with the one constructed from the 00-component of the energy- 
momentum tensor. In the spinorial representation this tensor has the form: 

= 2eohB9^^^9‘^^''^AB. (16) 

where eo is the permittivity of free space. Thus, the following two expressions must 
be equal: 


t 


fuv cl(k)cx{k) = 2eo / d-^r(pj^g9 9 <pAB- 


(17) 


The left hand side with the use of m can be expressed in terms of spinors as follows: 


t 


He 

Huj c\(k)c\{k) = — 




(18) 


Comparing this result with (HID we must choose 'Ye = \/hc(2eo and we obtain finally: 

d^k 


(pABix) = \/ ?ic/2eo 




kabb 


c+(fc)e-*'='“ -f cL(fc)e*'='^ 


(19) 
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In the standard approach to linearized gravity (see, for example, i) one starts 
from the decomposition of the metric tensor into the background metric (usually 
Minkowskian) and a small addition, g^i, = rj^i, + Note that the smallness of 
does not necessarily imply the smallness of its derivatives so that the linearized 
Riemann tensor does not have to be small. Next, one proceeds to express the linearized 
Riemann tensor (I shall keep referring to this tensor as the Riemann tensor even though 
in this case it effectively reduces to the Weyl tensor) in terms of and its derivatives. 
I will bypass all these intermediate steps and following [5] I will connect directly R^^Xp 
to its representation by the symmetric fourth-rank spinor, 

<fABCD = (20) 

i (sJ^S^;f^cj,ABCD + ■ (21) 


In order to express the field operator RpvXp in terms of annihilation and creation 
operators of gravitons we need the normalization factor yc in the formula: 


^abcd{x) = "fa 


dfk 


(27r)3/2fc 


kakbkci^d 


g+ik)e 


-ik-x 


9-ik)e 


ik-x 


( 22 ) 


The calculation of yc is a nontrivial task since the expression for the energy of the 
gravitational field in terms of the Riemann tensor is yet to be determined. Therefore, 
at this stage I will use a very good substitute — the Bel-Robinson tensor. The time 
component of this tensor is positive and it may play sometimes the role of 

the energy. It has even been called the superenergy [ini nn ng. As is shown in 
[Appendix B[ the calculation of yc in this way gives: 

yc = V^lp, (23) 

where (.p = sJhGjc? is the Planck length. Of course, from pure dimensional analysis 
we can deduce that yc has the dimension of length so that it must be proportional 
to the Planck length but the determination of the numerical coefficient requires the 
calculations in [Appendix B| Having determined yc we may freely pass from the 
description in terms of oscillators labeled by k and A to the description in terms 
fields living in space-time with the use of the formula: 


4>abcd{x) = x/SttHCJc^ 


(Pk 


(27r)^/^/c 


kakbkckd 


9+ik)' 


-ik-: 


9-{k)t 


(24) 


5. Wigner functional of the electromagnetic field 


The main tool in the study of the fluctuations will be here the Wigner functional 
at finite temperature. This functional for the electromagnetic field was obtained in 
[3] but I shall sketch here the main steps of the derivation to continue easily to the 
gravitational case. I start with the Wigner function at finite temperature for the 
one-dimensional harmonic oscillator [T3l[l4], 


Wt{xtp) = Cexp 


—2 tanh 


hu! 

2 W 


H{p, x) 
tujj 


(25) 


where ff(p, x) = p^/2m-\-muj‘^x^/2 is the Hamiltonian of the harmonic oscillator. The 
ratio N = H{p,x)/fkjj is the number of quanta (the energy divided by the energy of 
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one quantum) expressed in terms of classical variables (p, x). At T = 0, i.e. in the 
ground state, Wq = C exp{—2N). The normalization constant C is unimportant since 
for the infinite number of oscillators only the relative probabilities can be determined. 

The Wigner functional of the thermal state of the electromagnetic field is 
constructed by replacing a single oscillator with the whole collection of oscillators 
labeled by k and A [3], 




— j— r'\) (^E{r )-E{r') + (?B{r )-B{r')) 


■ (26) 


The electromagnetic correlation function /^(r) is the three-dimensional Fourier 
transform of the function appearing in the Wigner function for the one-dimensional 
oscillator, 




(flk tanh(£Qfc/ 2 ) 
(27r)3 k 


ihr 


eo 


2'KhciQr su\U{'Kr/(. q) ’ 


(27) 


where the quantum thermal length is (.q = ?ic//cBT=0.0023m/T[K]. In the derivation 
of (1^ I used (EH) and the following relation: 


= - {E{r)-E{r') + c^B{r)-B{r')) . (28) 

From the integral representation (l?7l) of /e(x) we obtain in the classical limit, 

lim/ e(x) = ^(5(r). (29) 

fi—>-0 z 

Therefore, in this limit the Wigner functional becomes the Boltzmann distribution, 
WiM = eM-EEM/kBT). (30) 

At the other end, at T = 0, the Wigner functional is equal to: 


^EM 


E, B] = exp 


-eo d-^r 


^^^,E{r)-E{r')+^B{r)-B{r') 


(31) 


2'K'^hc\'i — r 'P 

i.e. equal to exp(—2iV) where N is the total number of photons as calculated by 
Zeldovich [II1I3]- 


6. Wigner functional of the gravitational field 


The thermal Wigner functional of the gravitational field can be constructed in the 
same way as for the electromagnetic field. The resulting formula is: 


Wq (i?) = exp 




(32) 


where £ij = Riojo and Bij = ^CikiRjQ^’’ are the so called electric and magnetic parts 
of the curvature tensor [Hill]. In the derivation of this formula I used the following 
gravitational counterpart of (l28l) derived with the help of EH: 

4’ABCDi'<^)9°^'^9°^^9°^^9°°°(l}ABCDir') = (^*j (»’)^y (’’') + . (33) 

ij 

The gravitational correlation function fair) can be explicitly evaluated and is given 
by the following counterpart of the formula (071) . 


Mr) = T 


d^k 2tanh(^Q/c/2)^jy;,.^ 


(27r) 


kf' 


1 


STT^fgr 


y -f ln(C) ln(l + C) + Li 2 (l - C) + Li2(-0 


(34) 
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where (.q = GhsT/c'^ = 1.14 x 10“®^mT[K] is the gravitational thermal length, 
f = coth(7rr/2t!Q), and Li 2 is the dilogarithm function. 

There is a great similarity between the probability distributions of various field 
configurations for the electromagnetic field and the gravitational field. In both 
cases the smaller the distance between the points, the more likely it is that the 
electromagnetic field or the curvature tensor at these points will have opposite signs. 
Thus, the formula (l32ll may be viewed as a realization of the Wheeler concept of the 
virtual gravitational foam [18] . Of course, in the linearized version of gravity there is 
no room for fluctuations of topology envisaged by Wheeler but the rapid changes of 
geometry at small distances do show a certain similarity with Wheeler’s foam. 

In the classical limit in the electromagnetic case we obtained the standard 
Boltzmann distribution dsni) and the same result holds in the gravitational case, 
Wq{R) = exp{—EG/kBT), where 


Eg 



£ij jr )£ij (r ') + % (r {r ') 

32Tr^G\r — r'\ 


(35) 


One may check that Eg is indeed equal to o by expressing dSSl), with the use of dm 
and (1331) . in terms of annihilation and creation operators. 

Eg =^hujgl{k)gxik). (36) 

The equality of these two expressions proves the self-consistency of the whole 
procedure. The classical (high-temperature) limit gives the correct expression for the 
energy both in the electromagnetic and in the gravitational case. The nonlocal form of 
Eg fully confirms the belief that there is no “local gravitational energy-momentum” 

0 - 

There are, however, substantial differences between the gravitational and electro¬ 
magnetic correlation functions. At large distances the gravitational correlation 
function /G(r) does not fall-off exponentially, as in electromagnetism, but has a long 
tail equal to its classical limit /g(?') ~ /g(’') ~ l/(327r^£G’")- There is also a striking 
difference in the behavior at T —?> 0, as described in the next Section. 


7. Problems with the ground state 

A puzzling phenomenon is the logarithmic divergence of the gravitational correlation 
function /G(r) in the limit, when T —^ 0. This fact indicates that there might be a 
problem with the gravitational ground state. 

The ground state of the quantized gravitational field has been studied before by 
Kuchaf m and Hartle [5D] . They arrived at a formula for the ground state functional 
but their analogy between the electromagnetism and gravity is based on a different 
paradigm. They saw the analogy between the electromagnetic potentials and the 
metric tensor while in my approach, based on the spinorial description, there is a 
correspondence between the electromagnetic field tensor and the Riemann tensor. 

The roots of the problem are in the expression (13511 for the energy. To see 
the connection between the energy and the ground state let us consider the wave 
functions of a one-dimensional harmonic oscillator in the position representation and 
the momentum representation: 
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In the position representation the exponent is equal to the potential energy divided 
by hjjj and in the momentum representation it is the kinetic energy divided by huj. 
Applying the same prescription to the electromagnetic field we obtain the Wheeler 
expression [TH] (also rederived in [TH] and [10]). This follows directly from the 
observation that the division by hio in quantum mechanics is represented in quantum 
field theory by the following kernel 


J huk'' hc\r — r'p ’ 

However, in the gravitational case we do not reproduce the Kuchaf-Hartle result but 
we encounter the same logarithmic divergence that appeared in the Wigner functional. 
The convolution of this kernel with the already present kernel l/\r — r'\ in the formula 
(II3 for the energy leads to the logarithmic divergence at large values of r. The same 
divergence appears already in the formula for the number of gravitons obtained by 
setting h{k) = 1 and H — 2 in (TTII) . 

The divergence of the gravitational probability at T = 0 can also be related to the 
behavior of the two-point function (sometimes called Riemann correlator [21j l at small 
distances. With the use of the formulas (nni) and (EB) we can evaluate the vacuum 
expectation value 


mAB...Lix)hB...Li^'m (39) 

for the spinor fields of rank 2H. The resulting expression will have 2H derivatives 
acting on A(+^(a; — x') —the invariant (positive frequency) solution of the wave 
equation. For a Gaussian distribution there is a simple relation between the matrix of 
the second moments (xiXj) and the matrix in the probability distribution. Namely, one 
is the inverse of the other. There exists an analogy between this discrete case and our 
field-theoretic continuous case. The two-point function is the analog of the covariance 
matrix and the correlation functions Je and fc in (|Mll and (l32|) are the analogs of 
the inverse of the covariance matrix. The divergence of fc at zero temperature can 
be traced to the nonexistence of the inverse of the gravitational two-point function. 
The presence of four derivatives in this function results in the fourth power of k in the 
Fourier transform and its inverse is not integrable at small k. In the electromagnetic 
case we have the second derivatives and this leads to an integrable Fourier transform. 

A technical cure for the problem with non-integrability could be to disallow the 
interchange of the integrations over r and r' in the formula (l32)l and the evaluation of 
the limit when T —>■ 0. Unfortunately, the limit would then depend on the behavior of 
the Riemann tensor at large distances. Since the divergence is of the “infrared type” 
it can presumably be cured by assuming a finite radius of the Universe but this is 
outside the scope of this paper. The problem with the ground state might also mean, 
however, that there is some truth in Dyson’s hypothesis. 


8. Conclusions 

I have shown that the quantization of the linearized gravitational field that employs 
only the Riemann tensor (and not the metric tensor) can be achieved without any 
reference to the canonical formalism. In this approach the complications arising in 
the process of extracting true degrees of freedom never appear. The main result is the 
evaluation of the probability distribution of various geometries in the thermal state 
of the gravitational field. An encountered puzzle is the divergence of this probability 
distribution at T —>■ 0 that might signify some problems with the ground state. 
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Appendix A. Proof of the equality (jllll 


This equality can be proven by first inverting the Fourier transform in (11011 at t = 0 
which leads to: 


ka^b ■ ■ ■ KLa+{k) + kakb ■ ■ ■ kLal (— fc ) = - 


tUr 


e-^^-^<fAB...L{r), (A.l) 


7 J (27r)3/2 

where KA{k) = KA{—k). In the next step I evaluate the product of both sides with 
their hermitian conjugate inserting the unit matrices to enable the summation 

convention by balancing dotted and undotted indices. Using the relations: 

(A.2) 


.,0.4.4^ 


.,0 AA: 


k-a 9''^‘^ KA = k = k^g^^^^^KA, 


^OAA? 


.,0.4.4^ 


Ka = 0 = k^g^^‘^‘KA, 


we eliminate the cross terms and obtain: 


a\_{k)a+{k) + al(—fe)a_(—fc) = 


1 


j2k2H-2 (27r)3 


d'^r e 


ikr 


d^r' 


X 0)g°^^g°^^... g°^^^AB...L{r', 0). (A.3) 

Multiplying both sides of this equality by h{k) and integrating over k we obtain the 
relation (El- 


Appendix B. Bel-Robinson tensor 

To find the normalization factor jg I use the Bel-Robinson tensor [giiniEiin]. tms 
is a fourth-rank tensor which to some extent can play the role of the energy- 

momentum tensor; in the Einstein-Maxwell system it satisfies the continuity equation 
= 0. Since the sum of the contributions from electromagnetism 
and gravity is conserved, the normalization of the electromagnetic part determines the 
normalization of the gravitational part. To evaluate 7 g I shall use the spinorial form 
of the Bel-Robinson tensor [6] , 

rppi^Xp _ ^0 nAAuBB \cc nbo 

^EM ~ 9 y y y 

^ „4 .... 

rrPvXp _ „iJ.AA ^vBB „XCC „pDD j, j, /td 

9 g ~ IQ-j^q 9 y y y VaBCD'PaBCD, (R-Jj 

where S means the symmetrization with respect to the listed indices. The electro¬ 
magnetic and the gravitational contributions to the super-energy can be evaluated 
with the use of an and dm, 

Jd\nT = l'^f^k^ci{k)cx{k), (B.3) 

Jd\T^°°° = J^^hu;k^gl{k)g^{k). (B.4) 

These contributions must have the same form since they both represent the same 
physical quantity. Therefore, 7 g = \/^fp, where £p = •JhGjl? is the Planck length. 
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